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Abstract: Zlatar Ski Center has an extremely high potential as a ski resort that could be 
among the most visited ones in Serbia. What characterizes this center is the proximity of 
Kokin Brod Lake, the Tara River and the like. Therefore, Zlatar should beconsidered both as 
a ski center and in terms of other diversities. The aim of this paper is to analyze the 
possibility of expanding thecapacity of skiresorts by increasing the possible flow of skiers 
on new ski lifts and ski slopes. The methods used in this paper belong to operational 
research, specifically, the geometric interpretation of linear programming. It is also shown 
howto optimize the extended capacities. If the capacities of skislopes are to be increased, it 
would initiate greater investment in accommodation capacities in the area of the municipality 
of Nova Varo, which would affect its improvement in terms of tourism. 
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Primena linearnog programiranja na dimenzioniranje 
proSirenog kapaciteta ski centra Zlatar u Novoj 
VaroSi 


SaZetak: Ski centar Zlatar poseduje izuzetno vis ok potenciyjal kao skyaliste koje bi moglo da 
bude u rangu najposeéenyih u Srbiji. Ono Sto ovaj centar isti¢e, 1 u pozitivnom smislu 
izdvaja od drugih, je blizina jezera Kokin Brod, reke Tare i sli¢no. Prema tome, Zlatar ne bi 
trebalo razmatratisamo kao skicentar, vec isa aspekta posedovanja i drugaciyih diverziteta. 
Cil rada je analiza mogucnostiproSirenja kapaciteta skyali8ta povecanjemmogu¢éeg protoka 
skijaSa na zi¢arama i kapaciteta zi¢ara. U tom smislu, predvida se i povecanje broja staza, 
kao 1 interkonekcya izmedu njih. Metode koje su kori8¢ene u radu pripadaju operacionim 
istrazivanjima, konkretno, geometrijs ka interpretacija linearnog programiranja. Prikazano je 
da je na ovaj nacin moguée optimizovati proSirene kapacitete. Povecéani kapacitet ski staza 
inicirao bi vece ulaganje u sme$ tajnekapacitetena prostoru opStine Nova Varo§, cime bise 
uticalo na njeno unapredenje u turisti¢kom smislu. 
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Kljuéne reéi: lnearno programiranje, ski liftovi i ski staze, optimizacija proSirenog 
kapaciteta 
JEL klasifikacija: C44, Z32 


1. Introduction 


The paper provides an approach to establishing the required number of lifts in the case of 
expansion of Zlatar Ski Center assuming an increase in the number of skislopes. The starting 
point is that the chair lift that supplies more skislopes must not haveacapacity greater than 
the capacity of the skislope with the lowest capacity that gravitates to that lift. The reason is 
a security aspect (Wheeler, 2012). In the event that some of the skislopes are closed and the 
capacity of the chair lift to which these ski slopes gravitate is fulfilled, skiers would be 
directed to one or more of the remaining ski slopes in service. In that case, the open ski 
slopes that are in service should not beoverloaded. The flow of skiers should not be higher 
than the maximum availability of the ski slopes (Poulhés & Miria, 2017; Revilloud et al., 
2011). For this reason, the sizing of the chair lifts is realized in the stated way. In the 
following chapters a model for sizing the capacity of lifts in accordance with the existing 
system of ski slopes and the projected capacity of new ones (Wheeler, 2012) will be 
provided. 


The main hypothesis of this research is that it is possible to resize the capacity of Zlatar Ski 
Center, and based on mathematical model and linear programming implementation the 
optimal solution of extended model dimensioning of ski slopes and ski lifts could be 
calculated. 


Optimization of tourist capacities with the necessary logistics is permanently present topic 
(Jangra & Kaushik, 2021), and the application of the linear programming for that purpose 
was considered in papers published before (Borrelli et al., 2003; Kniff & Oosterhaven, 
1990). 


2. Mathematical models for ski slopes capacity dimensioning 


The capacity ofaskilifts systemis projected by applying flow analysis and Kirchhoff rule 
that total input to the systemis equal to the output fromit. It is necessary to respect the 
equation of flow continuity (Madarevic, 1969). The flow of skiers on the ski lifts is the 
function of the capacity of skislopes by whichthe systemis getting relieved. Chair ski lifts 
that lead to the peaks havethe largest capacity. They are directed from the accommodation 
facilities, like hotel complex or ski bus terminus to the highest peaks. From these points, 
there is aflow divergence to smaller capacity T — bar type skilifts, which cover the complex 
of skislopes (Markovié, 2019a; Suvajdzi¢, 1972). According tothe fact that the organization 
of skiers’ flow is the main activity in the considered area (Zrni¢, 1993), ski center represents 
a big system of internal transport. 


2.1. Chair lift dimensioning 


The capacity of ski lifts is dimensioned on the basis of the required flow of skiers, i.e. the 
available accommodation capacity of the hotel or the capacity of the bus line by which the 
skiers are brought to the starting point of the chair lifts. The capacity of the ski lifts leading 
from the starting point is the largest and matches the sumof capacities of smaller ski lifts 
(often T-bar type) which gravitate to the ski slopes capacity. According to that, the 
possibility of shifting skiers fromone skislope to the other has also to be considered. It is 
important that percentage of such flows must not exceed the maximum allowed ski slope 
capacity. Presentation of sucha situation is provided fora hy pothetical ski center in Table 1 
and Figure | (Markovic, 2019a). 
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Table 1 represents the schema for the share of skiers per hour crossing fromone ski slope 
(rows) to the other (columns). Values in the cells are expressed in percentage of the 
maximum capacity of the skislope, in skiers per hour that are expected to cross from one 
slope to another. 


Table 1: Share of skiers perhour(%) crossings fromone skislope to another and ski slopes 
capacity of a hypothetic ski center 


Capacity 
Skislops | 1 | 2 | 3 | 4 5 | 6/7 | 8 [skiers/h] 





eS a 


Source: Author’s research 


For example, the estimated share percentage of skiers crossing from the first ski slope to 
Slope 4, 5 and 7 is maximum 10%, while only 5% could be transferred to Slope 7. 
Accordingly, the capacity of each skis lope on which crossing is assumed has to be boosted 
according to enlarged number of skiers per hour. A similar approach is to be applied for 
other skislopes (Markovic, 2019a). In the case of designing a new ski slope, No. 8, which 
gravitates to the chair lift 2 from direction | and to the chair lift 1 fromthe same direction, 
including the inflow of skiers from other ski slopes, the layout will be as presented by 
scheme, Figure 1 (Markovié, 2019a). 


Figure 1: Ski slope 8 loading, by chair lifts and by adjacent ski slopes 
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Source: Author’s research 
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Required capacity on the jof the total: n directions fromthe center of the ski resort, n = 3; L, 
II and III direction, to the peaks is as follows: 


Q, [skiers/ h]< Ya; 


i=l 


? 


where: 
- gj [Skiers/h] — capacity of the 7 ski lift in the / direction, 
- m-— total number of ski lifts in the 7 direction. 


The total flow of all ski lifts in skiers per hour is expressed by: 


n m 


DZ, = 4 ; 


j=l i=l 
2.2. Ski slopes dimensioning 


Necessary capacity t, of the skislope g represents the projected hourly flow of skiers on that 
ski slope in the form: 


I S 
be [skiers /h] = D2. [skiers /h] + 2 Pe -Z; [skiers / h] 
72 am 
where: 
- s—totalnumber of skislopes, 
-  Z, [skiers/h]— projected hourly flow of skiers on ski lifts to which the g skislope 
gravitates, 
- — p,[%/100] — the share of capacity expressed through the hourly flow of skiers 
crossing fromi ski slope to the g skislope, 
- Z,[skiers/h]— projected hourly flow of skiers on i skislope, 
- _-totalnumber of skilifts by which the skislope gravitates. 


The capacity of all ski slopes is in the form: 
totalNo of all ski slopes 
T|skiers /h|= pa: 


g=l 
3. Design of necessary chairlifts and ski slopes capacity 


Since the aim of this paper is analysis of possible capacity expansion of Zlatar Ski Center, it 
is necessary to presentthe current situation first. Zlatar Ski Center current state is represented 
in Figure 2 with ski slopes I and II, as wellas the four-seat chairlift that leads from foot to 
the top, together with T—bar ski lift that leads from the fork point b, to the peak point a). 
Peak point a, is top point to which it is possible to get fromthe foot of the mountain as well 
as from the middle fork point: b;, Figure 2. This paper estimates that in addition to the 
mentioned peak point and middle fork point b;, new peak points should be designed: a> and 
a3, Figure 2. Skiers will be able to get to these points fromthe new base point S by newly 
designed chairlifts: 2,3 and 4 as suggested in Figure 2. Skiers are expected to access several 
newly designed skislopes, i.e. III to VI from top points a, and az, and fromthose new slopes 
middle points b> to b;, Figure 2, skiers would be able to get to the top again by newly 
designed smaller capacity T-bar type ski lifts. According to that, skiers are not forced to 
descend to the base point in orderto reach the top points again by high capacity ski lifts, but 
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with smaller capacity, T—bar type ones. Considering the unused benefits of this ski center 
due to high snow amount and average number of snow days during the year, Table 2, this 
paper proposes a solution tothe possible capacity expansion of this skicenter by introducing 
new fork points, new ski slopes and ski lifts as described. 


Table 2: Average snowfall at Zlatar 


Snow amount (week) Snow days (week) 
2.8 days 


3.3 days 
35 days 


March 23 cm 3.4 days 
1.6 days 


Source: Snow-forecast (2021) 





According to the abovementioned, newly designed T—bar ski lifts would be used for taking 
skiers from mid-points b;,i=2,... n,n =5 to pointsa;,i= 2, 3 on the peakas they will not be 
forced to take chairlift from the starting point S, Figure 2. 


Altogether, the existing and newly designed layout of the Zlatar Ski Center is shown in 
Figure 2. The full lines show skislopes, the dashed lines represent newly designed chairlifts 
in the directions toward a, anda; points. The Arabic numerals refer to chairlifts while Rome 
numerals refer to ski slopes. 


Figure 2: Zlatar Ski Center, existing and newly designed layout 





Legend 


I, Il... - skislope sign, 
1,2... - the sign of the 
chairlift to the first 
forks - a), a2, and a3, 
—— - existed skislope, 
== =- chairlift, new, 
== - ski crossway, 
== - existing chair lift 
=—=- ski crossway to 
fork point b,, 
- crossway from 
- rise chairlift to slope, 
La | bs} by-=bs ==". newly designed 
New fork points b;, ski slope, 
- “ski slope which is 
yet to be designed, 
=== - road to ski bus, 
" §  - newly designed 
start point. 


In order to increase the number of ski slopes and the chair lifts of Zlatar Ski Center, a new 
starting pointis planned towards the fork points ay and a3, as well as the crossing from one 
skislope to another, Table 1 and Figure 2. 
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Newly designed concept introduces the expansion of the capacity of the Center, where it is 
necessary to harmonize the capacity of the newly designed ski slopes and the required 
capacity of the newly designed chair skilifts, which leads to the optimization of the existing 
conditions. The required capacity of each of the fournewly designed chair lifts that supply 
five ski slopes of the newly designed Ski center Zlatar is established. Ski center is to be 
designed with the ski slopes capacity of 1000, 2000, 3000, 2000 and 4000 skiers per hour. 


The criterion for determining the required capacity of the ski lift is the fulfillment of the 
conditions that the capacity of the ski lift must be at maximum, but not higher than the 
capacity of the skislope with the lowestthrough pass, comparing the through passes of all 
skislopes that gravitate to that ski lift. 


“ce 


On the other hand, the capacities of chair lifts that supply a particular peak point “a”, Figure 
2, must be approximately equal to the capacity of all smaller capacity T-bar ski lifts that 
gravitate to the point “a”. Thecanacity of chair lifts is in correlation to the canacity of the ski 
slopes that gravitate to them and is to be established in chapters that follow. 


The required capacity of each of the newly designed chair lifts, including theexisting one, is 
being established. The calculation is also realized in compliance with the principle that the 
total capacity of chair lifts that supply skislopes mustnotbe greater than the capacity of the 
lowest capacity skis lope fromthe aspect of the hourly flow of skiers. When performing the 
calculation (Markovic, 2019a), the shifting of skiers fromone ski slope to another is not 
foreseen while the chair lifts supply the ski slopes in the following way: ski lift No. 1 
supplies the skislopes | and 2; skilift No. 2, supplies the skislopes 2,4, and 5; ski lift No. 3, 
supplies the ski slopes 3, 4 and 5; and ski lift No. 4, supplies the ski slope 5 as well as the 
newly designed one, i.e. No. 6, Figure 3. 


Figure 3: Ski center scheme 





Ski 
slopes: 
n=6 


Charrlifts 
m=4 


Source: Author’s research 


It is not allowed to achieve influxof skiers on the skislope greater than its capacity. On the 
basis of this, through the geometric interpretation of linear programming (Nadler, 1973; 
Markovic, 2019b; Borrelli et al., 2003; Kasana & Kumar, 2004), inequalities are formed for 
the maximum capacity (Nash & Sofer, 1996; Sarbjit, 2018) ofeach skislope:i, i =1, ... ,n=5, 
in the form: 


C, < 1000, for ski slope 1, 
C, + C, $2000 => C, < 2000 — C,, for ski slope 2, (1) 
C3 < 3000, for ski slope 3, 
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C,+ C3 < 2000, for ski slope 4, 
> C,< 2000—-C,,C,=0 > 2000 -C,>0>C, < 2000, (2) 


C,+C34+C, < 4000, for ski slope 5, 
C, < 4000 —(C,+C3),C, >0 => 
4000 -—C,—-C,20,C,2=0>C, < 4000. (3) 


Inequalities follow in the next form: 


C,<1000, C,<1000, and C;> 0, i=1, ...,m,m=4. (4) 


In this case C;,i= 1, ..., m, m =4, representthe capacities of the ski lifts (Dedijer, 1983) that 
should weigh the highest possible value, whereby it is not possible to supply the ski slopes 
with the flow of skiers that exceed the capacity of the skislope. It has to be noted that the 
requirement is that the capacity of the ski slope must not be exceeded even in the case of 
closing one of the ski slopes that are supplied by the same chair lift. 


The objective function that represents the maximum flow of skiers on allski slopes of the ski 
resort, and should strive for the maximum value is in the form: 


n=4 
0<F=)C,i=1+4,F, =C,+C,+C,+C,>max | (5) 
i=l 


Hence: 





F,=C,+C, +(2000-C,)+(4000-(C, +C,))= 
=C,+C, +(2000-C,)+(4000-(C, + 2000-C, )) 





F. =C,+C, +(2000-C, )+ 2000 = C, +4000, F. > 0. 6) 


Inequalities (1), (2), (3) and (4) are entered in diagram C,; — C2, Figure 4. A vertex 
coordinates of M,, My, M,and M, represent capacities of chair lifts 1, 2, 3, and 4, Figure 4. 


Figure 4: Chair lifts capacity domain and capacity of objective function (6) 
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maxim is 







1p Co, C3, Cy >0 
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-4000 


Source: Author’s research 
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These vertexes form the closed polygonal line that bounds the domain for the chair lifts 
capacities in skiers per hour, ie. C1, C2, C3, C4, (Markovic, 2019a), and are as follows: 


M (Cis 


CrCl =1,2,3,4. 


M,(0, 0, 2000, 2000), 7,(0, 2000, 0, 2000), 
M,,(1000, 1000, 1000, 2000), 17, (1000, 0, 2000, 2000). 


(7) 


The largest traffic on the chair skilifts is realized in the vertexes M3 and M, (7). As it is not 
possible for traffic on the chair lift to be zero, by the reason of its abolition in that case, the 
values of coordinates in vertex M, are to be adopted as adequate. That means that the 
capacity of the chair lifts will be C,=1000, C,=1000, C3=1000 and C,=2000 skiers per hour 
(7). This state is to be considered as optimal. 


The coefficient of skislope utilization from the aspect of cableway capacity is in the form: 














pe OY 50k. Se oem: 
* "5 1000 2s 2000 
ea Eg SGU Bo a oo, 8) 
"5. 3000 «5 2000 
k= Sos a OT og 100%: 
" 5, 4000 


Where: S;[skiers/h], i=1, ..., 5 — capacity of i - ski slope. 


Authors of the paper created and ran the program SKIERS 6, based on Williams (2003), 
taking into account calculations presented in this chapter, (1) to (8). The res ults are presented 
in the chapter that follows. The difference is that in new cases the percentage of skiers 
crossing fromone skislope tothe other is also included in account scheme. This option is 
described in the hypothetical example in chapter 2. 


4. Results ofthe analysis ofski flows at Zlatar Ski Center, Nova Varo§ 


By applying the SKIERS 6 programpackage, the results listed in the files were obtained as 
follows: 


Total capacity of all ski slopes: 12000 [skiers/h] 

Total capacity of all chair ski lifts: 5000 [skiers/h] 

The percentage utilization of the ski slopes is 41.66666 [%] 

Minimum of the capacity necessary forthe skislope 1, calculated with no crossing to this ski 
slope from the others is 1000 [skiers/h] 

Minimum of the capacity necessary for the skislope 2, calculated with no crossing to this ski 
slope from the others is 2000 [skiers/h] 

Minimum of the capacity necessary for the skislope 3, calculated with no crossing to this ski 
slope from the others is 1000 [skiers/h] 

Minimum of the capacity necessary for the skislope 4, calculated with no crossing to this ski 
slope from the others is 2000 [skiers/h] 

Minimum of the capacity necessary forthe skislope 5, calculated with no crossing to this ski 
slope from the others is 4000 [skiers/h] 


Minimum of the capacity necessary forthe skislope 6, calculated with no crossing to this ski 
slope from the others is 2000 [skiers/h] 
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The scheme of the anticipated percentage of crossings fromeach skislope to all the others 
[skiers/h] is as follows: 


Ski slope No. 1 No. 2 No. 3 No.4 No.5 No. 6 
Ski slope No. 1 0 200 0 0 0 0 
Ski slope No. 2 0 0 0 0 0 0 
Ski slope No.3 0 0 0 0 0 0 
Ski slope No. 4 0 0 0 0 200 0 
Ski slope No.5 0 0 800 400 0 0 
Ski slope No. 6 0 0 0 400 0 0 


Minimum of the capacity necessary for the ski slope | is: 1200 [skiers/h] 
Efficiency coefficient of the ski slope No. 1 is: 7(1)= 1.2 

Minimum of the capacity necessary for the ski slope 2 is: 2000 [skiers /h] 
Efficiency coefficient of the ski slope No. 2 is: (2)= 1 

Minimum of the capacity necessary for the ski slope 3 is: 1000 [skiers /h] 
Efficiency coefficient of the ski slope No. 3 is: 7(3)= 0.3333333 

Minimum of the capacity necessary for the ski slope 4 is: 2200 [skiers /h] 
Efficiency coefficient of the ski slope No. 4 is: 4(4)= 1.1 

Minimum of the capacity necessary for the ski slope 5 is: 5200 [skiers /h] 
Efficiency coefficient of the ski slope No. 5 is: n(5)= 1.3 

Minimum of the capacity necessary for the ski slope 6 is: 2400 [skiers /h] 
Efficiency coefficient of the ski slope No. 6 is: n(6) = 1 

Minimum of the capacity necessary for the all ski slopes is: 14000 [skiers /h] 








Efficiency coefficient of all ski slopes fromthe aspect of the minimum necessary capacity is 
in the sum of 0.53. 


It is stated thaton the skislopes | and 5, the utilization coefficientis higher than one due to 
the anticipated possibility of skiers crossing fromone ski slope to another. This causes the 
necessary increase in the capacity of the ski slopes in order to reduce the utilization 
coefficients of the skislopes 1 and 5. The improved solution is presented in the scheme of 
the ski resort, Figure 5. 


Figure 5: Ski slopes capacities 


2400 |Ski slopes 


n=6 





Charr lifts 
m=4 


Source: Author’s research 


Based on Figure 5, through the geometric interpretation of linear programming, inequalities 
are formed for the maximum capacity of each ski slope: i, i =1, ... , n=6, in the form that 
follows: 


C, < 1200, for ski slope 1, 
C,+C, < 2000,C, =>0 => C, = 2000 — C, > C, < 2000, for ski slope 2, (9) 

C3 < 1000, for ski slope 3, 
108, 
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C,4+ C3 < 2200, for ski slope 4, 
> C3 < 2200—C,,C;=>0 > 2200 -—C, >0>C, < 2200, (10) 


C,+C34+ C, < 5200, for ski slope 5, 
C, < 5200 — (C,+C3),C, = 0 > 5200—C,—C,20,C,;>05C,—5200>05 
C, < 5200, 


C, < 2400, for ski slope 6. (11) 
According to inequalities (9), (10) and (11), there follows: 
C,<1200, Cr<2200, and C;= 0, i=l, ..., m= 4. (12) 


Here, C,, i=l, .... m=4, according to (12), representthe capacity of the chair lifts that should 
strive forthe highest possible value, while it is not possible to supply the ski slopes with a 
flow of skiers that exceeds the capacity of them. It should be noticed thatthe capacity of the 
ski slope should not be exceeded even in the case of closing one of the skislopes that is 
supplied by the same chair lift. 


The objective function thatrepresents the maximum flow of skiers on allski slopes of the ski 
resort, and should strive for the maximum value, is in the form: 


n=4 
O< Be Ca =1+4,F. =C,+C,+C,+C, > max , Figure 6. (13) 


i=l 


Figure 6: Chair ski lifts capacity domain 


The ¢ondition that has to 
be met is: 
1, Co, C3, Cy >0 





-5200 0 00 2000 C 
Source: Author’s research 


Vertex coordinates M,, M>, M3, M, of the domain, formed by closed polygonal line bounding 
the domain area that represents the chair lift capacities C,, Cs, C3, C4 [skiers/h], (12), are as 
follows: 


M, (Cy,Cigs Cas Cig) 8 =1,2,3, 4. 


il? 


M, (0, 0, 1000, 2400), 44,(0, 2000, 2000, 2400), (14) 
M,(1200, 800, 1000, 2400), 7, (1200, 0, 1000, 2400). 
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Objective function values in vertexes, according to (13), are as follows in (15): 


Fy, = 3400,F,y, = 4600,F.4, = 5400, F.y, = 4600. (15) 


5. Discussion 


According to previously mentioned in equation (14), the biggest flow on the chair lifts is 


assumedin vertex M3. That concerns the chair lifts capacities in the sizes: C\=1200, C,=800, 
C3=1000 and C,=2400 [skiers/h], (14). This state is considered to be the optimal one. 


Efficiency coefficients of the ski slopes according to chair lift capacities are as follows in 
(16): 


C, 1200 _ _C+C, _ 2000 

















was 1> 100%, k,,. =1> 100%, 
S, 1200 S, 2000 
oo i 0t te 2 OO pis sie, (16) 
S, 1000 S, 2200 
bg oe IT oo sien k= = a 00, 
| ca 5200 "65 2400 


Where: S; [skiers/h], i=1, ..., 5 — capacity of i- ski slope. 
6. Conclusion 


This paper presents a mathematical model, geometric interpretation of linear programming, 
as well as the Standard Maximum Method andits usage in the field of transport. The aim of 
the paper is to analyze the possibility of expanding the potential of Zlatar Ski Center by 
designing newskislopes as well as installing additional chair lifts in order to complete the 
tourist offer of that Center. Additionally, the paper also provides a model for finding the 
optimal capacity of additionally designed ski slopes. The ski center is studied as a large 
systemof internal transport, considering that the analysis of the flow is the main activity. 


The obtained results showed that it is possible to confirm the hypothesis stated at the 


beginning of this paper and get the optimal solution of the expanded capacity of the ski 
center by applying linear programming and the presented methodology. 


This paper alsopresents the application of a software program, developedby the authors, to 
establish the coefficient of utilization of the ski slopes as well as the entire s ki resort. The 
percentage of skiers who cross fromone skislopeto another is also covered, together with 
“goal” function ofthe problem and detecting its maximum. Finding the maximum capacity 
of the projected systemis also achieved. 
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